Prepare all of the needed data:

## Warning: Duplicated column names deduplicated: 'UTM' => 'UTM_1' [7]
## Parsed with column specification:
## cols(
##   .default = col_double(),
##   Group = col_character(),
##   GMU = col_character(),
##   Basin = col_character(),
##   Creek = col_character(),
##   UTM = col_character(),
##   Descriptor = col_character(),
##   `Hyb Notes` = col_character(),
##   `PVA Extinction` = col_character(),
##   `PVA Ext L95` = col_character(),
##   `PVA Ext U95` = col_character()
## )
## See spec(...) for full column specifications.
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular

Plots

Corrected vs. observed stats

## Warning: `group_by_()` is deprecated as of dplyr 0.7.0.
## Please use `group_by()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.

Theta trajectories